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First-Order  Parametric  Model  of  Reflectance  Spectra  for  Dyed  Fabrics 

Introduction 

Parametric  models  of  NIR  reflectivity  spectra  for  dyed  fabrics  are  presented,  which 
provides  for  both  their  inverse  and  direct  modeling1.  The  dyes  considered  contain  spectral 
features  that  are  of  interest  to  the  U.S.  Navy  for  the  purpose  of  simulating  NIR  spectra 
corresponding  to  various  dyes  and  mixtures  of  dyes  in  camouflage  textiles.  The  results  provide 
validation  of  the  constructed  parametric  models,  within  reasonable  error  tolerances,  and  will 
contribute  to  the  development  of  a  parameter  space  for  process  optimization  of  dyed  fabrics. 
Given  the  fabric  and  dyes  available,  an  optimization  algorithm  could  be  developed  that  predicts 
the  optimal  dying  process  and  fabric/dye  combinations  needed  to  provide  multi- spectral 
concealment  for  a  given  environment.  To  reach  the  point  of  process  optimization,  different 
absorbing  dyes,  fabrics,  and  parametric  reflectance  models  must  first  be  examined.  The  theory 
underlying  the  parametric  model,  as  well  as  the  dyes  and  fabric  chosen  for  this  analysis,  are 
described  in  subsequent  sections.  The  formulation  of  the  parametric  model  is  based  on  the 
inverse  analysis  approach1. 


Materials  and  Methods 

Dyes.  The  reflectance  characteristics  of  4  dyes  were  examined.  Each  dye  had  a  wavelength  of 
maximum  absorption  located  in  either  the  near  or  short-wave  infrared  bands,  with  appreciable 
spectral  features  extending  into  the  visible  spectrum2  (Table  1).  The  dyes  were  analyzed 
individually,  but  have  potential  to  be  used  in  combinations  to  create  custom  spectral  responses, 
i.e.  reflectance  or  transmission. 

Fabric.  The  base  fabric  considered  for  dye  application  was  a  50%  Nylon,  50%  Cotton  blend 
Ripstop  fabric,  i.e.  50/50  NYCO  Ripstop.  Ripstop  fabrics  are  woven  fabrics  with  reinforcement 
threads  interwoven  at  regular  intervals  in  a  crosshatch  pattern3.  The  fabric  was  characterized  as 
having  a  rough  surface,  finite  and  non-uniform  thickness,  heterogeneous  composition,  and  a  four 
color,  AOR1  desert  camouflage  pattern  design  (Figure  1). 

Dyed-Fabric  Sample  Preparation.  For  each  dye  and  corresponding  solvent  (Table  1),  a  stock 
solution  of  30  g  L"1  was  made.  Serial  dilutions  of  each  stock  solution  were  prepared,  consisting 
of  3.0,  1.5,  0.3,  0.15  and  0.03  g  L1.  Twenty  dye-solutions  were  prepared  in  total.  One  square 
yard  of  50/50  NYCO  Ripstop  fabric  was  cut  from  the  roll  and  placed  under  a  Freeman  Schwabe 
cutting  press.  A  2”x2”  Fremont  steel  rule  cutting  die  was  used  in  combination  with  the 
mechanical  press  to  make  20,  2”x2”  fabric  samples.  Dye-fabric  samples  were  prepared  by 
pipetting  1  mF  of  each  dye  solution  onto  the  center  of  each  2”x  2”  fabric  piece.  The  pipetting 
speed  was  adjusted  as  to  avoid  pooling  of  the  dye  solution  on  the  fabric  surface.  Dye-fabric 
samples  were  given  24  hours  to  dry  before  total  reflectance  measurements  were  made. 

Spectral  Measurements  of  Dyed-Fabric  Sample.  Measurements  of  total  reflectance  were  made 
using  a  UV/Vis/NIR/SWIR  spectrophotometer  (Perkin-Elmer®  Fambda™  1050,  Waltham,  MA, 
USA)  in  combination  with  UV  WinFab  (V6.0.3,  Perkin  Elmer)  software.  Reflectance  was 
measured  from  250  nm  to  2500  nm  at  5  nm  stepping  increments.  All  measurements  were  8 
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degree  hemispherical  in  total  reflectance  collection  mode.  The  Lambda  1050  incorporated  a 
double  beam,  150  mm  integrating  sphere  housing  a  photomultiplier  tube  (PMT)  detector  for  the 
UV-Vis  (175  -  860  nm)  region  and  an  Indium  Gallium  Arsenide  (InGaAs)  detector  for  the  NIR- 
SWIR  (860  -  2500  nm)  region.  Radiation  sources  included  a  deuterium  lamp  for  operation  in  the 
UV  (175  -  319  nm)  and  a  tungsten  halogen  lamp  for  use  in  the  Vis-NIR-SWIR  (319  -  3000  nm). 
A  light  source  change  at  319.2  nm  and  a  detector  change  at  860.8  nm  occurred  automatically 
during  monochromator  slewing.  For  UV-Vis  wavelengths,  the  slit  width  was  fixed  at  4  nm.  For 
NIR-SWIR  wavelengths,  the  slit  width  was  set  to  “Servo”  mode,  which  automatically  adjusted 
the  slit  width  during  scanning  to  maintain  constant  energy  at  the  detector. 

Parametric  Model  of  NIR  Reflectance  Spectra  for  Dyed  Fabrics 

Physically  consistent  parametric  models  provide  a  means  for  investigating  the  spectral 
response  of  a  given  system,  consisting  of  a  given  material  or  combination  of  materials,  and 
geometry,  in  terms  of  both  inverse  and  direct  analysis  [4, 5, 6, 7],  For  reflection  from  a  dyed  fabric 
of  finite  and  nonuniform  thickness,  whose  material  composition  is  heterogeneous,  a  general 
framework  for  parametric  modeling  of  reflectivity  spectra  can  be  in  terms  of  formulations  based 
on  the  Beer-Lambert  Law  or  its  approximation,  which  is  adopted  for  the  parametric  model 
considered  here  [8,9].  Accordingly,  a  parametric  model  based  on  characteristic  trend  features  of 
reflectance  spectra  from  fabrics,  which  assumes  a  first-order  absorption  process  within  an 
average  penetration  depth,  is 


where 


and 


^,m(A)  =  (1-/(A))R/(A) 
/(A)  =  2  A, a, (A) 


(1) 

(2) 

(3) 


where  a,  is  the  absorption  coefficient  of  component  j.  The  formal  procedure  underlying  inverse 
analysis  based  on  Eqs.(l)-(3)  entails  adjustment  of  the  absorption  coefficients  defined  over  the 
entire  wavelength  region,  if  assumed  unknown,  and  the  coefficients  Aj.  This  approach  defines  an 
optimization  procedure  where  the  reflectivity  spanning  the  range  of  wavelengths  is  adopted  as 
the  quantity  to  be  optimized.  Constraint  conditions  are  imposed  on  the  reflectivity  by  minimizing 
the  objective  functions  defined  by  Eq.(3),  respectively,  where  Rm(An)  is  the  measured  or  target 
reflectivity  for  wavelength  A„.  The  quantities  wn(n=l,...,N)  are  weight  coefficients  that  specify 
relative  levels  of  influence  associated  with  constraint  conditions  Rm(An). 

The  mathematical  foundation  of  inverse  analysis  based  on  Eqs.(l)-(3)  is  that  of  least- 
squares  parameter  optimization  [10,11].  A  significant  aspect  of  least-squares  parameter 
optimization  is  the  choice  of  a  sufficiently  complete  set  of  basis  functions.  This  implies  that  all 
possible  modes  of  a  given  process  can  be  modeled  parametrically  by  linear  combinations  of  these 
functions.  Accordingly,  the  parametric  models  defined  by  Eqs.  (l)-(3)  adopt  absorption 
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coefficients  as  basis  functions.  The  general  form  of  Eqs.  (1)  and  (2)  are  based  on  trend  features 
characteristic  of  reflectance  spectra.  Accordingly,  the  parametric  model  (Eqs.(l)-(2))  is  based  on 
a  first-order  extinction  model  of  reflectance,  whose  general  form  is  consistent  with  spectral 
features  of  diffuse  reflectance.  A  derivation  of  Eq.(l)  follows  from  the  first-order  relation  for 
extinction 


(4a) 

where 


R f+d 

Rf 


a  y.aeXM) 


aext  (h)-N  (C,  +  C,  ) 


(4b) 


and  N  ,  Cabs,  Cscat  are  the  particle  density,  absorption  and  scattering  cross  sections,  respectively, 
with  the  conditions  that 


and 


Rf,j 

Rd=-^ 

R, 


M 


(5) 

(6) 


where  Rf,  Rd  and  Rf+d  are  the  reflectivities  of  the  base  fabric,  pure  dye  and  dyed  fabric, 
respectively,  and  C;  and  aj  are  the  concentrations  and  absorption  coefficients  of  the  jth  dye. 
Equation  (5)  follows  from  the  assumption  that  the  contribution  to  the  reflectivity  due  to  the 
fabric  can  be  factored  out  from  that  of  the  combined  system  of  dye  and  fabric.  The 
proportionality  condition  given  by  Eq.(6)  represents  an  important  aspect  of  the  constructed 
parametric  model.  This  condition  is  necessary  in  that  the  extinction  coefficient  Eq.(4b)  is  based 
on  absorption  and  scattering  coefficients,  Cabs  and  Cscat,  respectively,  with  respect  to 
transmission  through  a  finite  thickness  of  particulate  material,  and  thus  does  not  provide  for  a 
well  defined  representation  of  extinction  resulting  from  both  transmission  and  reflection. 

Results  and  Discussion 

A  parametric  model  was  examined  to  characterize  the  diffuse  reflectance  properties  of 
dyed  fabrics.  The  model  was  constructed  using  a  first-order  relation  for  extinction  Eq.(4a). 
Consistent  with  this  relation,  diffusely  reflected  radiation  is  related  to  changes  in  scattering, 
absorption,  and  thickness  of  a  dyed  fabric  layer.  The  layer  is  assumed  infinitely  thick,  leaving 
reflectance  in  terms  of  the  relative  proportions  between  absorption  and  scattering.  Due  to  the 
complex  nature  of  scattering,  and  the  uncertainty  in  microstructure  and  dye  molecule  orientation, 
the  combined  influence  of  absorption  and  scattering  are  represented  phenomenologically,  and 
thus  not  associated  directly  with  any  physical  processes.  Accordingly,  the  function  f(X)  defined 
by  Eq.(2)  represents  the  combined  influence  of  all  underlying  physical  processes  contributing  to 
reflectance  within  the  dyed  fabric  layer,  i.e.,  reflection,  refraction,  diffraction,  absorption.  The 
function  f(X)  is  a  lumped  parametric  representation  consisting  of  a  known  absorption  coefficient, 
a,  intrinsic  to  each  dye  and  a  parameter  A.  The  absorption  coefficient  is  determined  through 
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inverse  analysis  of  transmission  spectra  for  the  dye  in  solution.  The  parameter  A  is  adjusted  with 
respect  to  Eq.(3),  i.e.,  objective  function  minimization,  and  accounts  for  statistical  variations  in 
fabric  thickness,  dye  distribution  in  fabric,  and  variations  in  concentration.  The  Rj(X)  denotes  the 
reflectance  of  the  untreated  fabric.  It  is  important  to  note  that  the  absorption  coefficient,  a,  is  a 
known  material  property  as  opposed  to  a  phenomenological  parameter,  i.e.,  A. 

Figures  (2)-(5)  show  reflectance  spectra  of  the  four  dyes  under  analysis  (Table  1)  as  a 
function  of  dye  concentration  and  wavelength.  The  reflectance  spectrum  of  a  control  fabric 
containing  no  dye  is  shown  in  each  figure  for  comparison.  Each  experimental  spectrum  was 
modeled  by  adjusting  parameter  A  until  the  objective  function  (Eq.  (3))  was  minimized.  Figures 
(6)-(9)  show  comparisons  between  modeled  and  experimental  data  for  each  dye  concentration. 
The  good  consistency  of  these  fits  demonstrates  the  acceptable  prediction  accuracy  of  the 
parametric  model  for  diffuse  reflectance  by  the  dyed  fabrics  considered.  The  results  also  show 
that  absorption  coefficients  for  dyes,  which  have  been  determined  by  inverse  analysis  of 
transmission  spectra  for  dyes  in  solvents,  can  be  validated  as  reasonable  estimates  of  absorption 
coefficients  for  dyes  in  fabrics.  References  [12,13]  describe  the  inverse  analysis  procedure  for 
determination  of  the  dye  absorption  coefficients  considered  in  this  study.  Supplementary 
information  concerning  the  estimated  dielectric  properties  of  these  dyes  as  determined  by  inverse 
analysis  is  given  in  the  appendix. 

Conclusion 

The  fitting  of  modeled  to  experimental  reflectance  spectra  shown  in  Figures  (6)-(9),  as  well 
as  the  relatively  good  consistency  of  these  fits,  establishes  proof  of  concept  for  constructing 
mappings  from  a  process-parameter  space,  i.e.,  dye  concentrations,  to  a  model-parameter  space, 

i.e.,  parameters  Aj  for  reflectance  spectra.  This  in  turn  establishes  that  the  parametric  model 
defined  by  Eqs.  (1)  and  (2)  can  be  adopted  for  prediction  of  dyed  fabric  reflectance 
characteristics  for  process  parameters  that  can  be  mapped  into  a  physically  consistent  range  of 
model  parameters. 
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Table  1:  NIR  and  SWIR  Absorbing  Dyes 


Dye  Name 

Absorption  Max  (  nm) 

Solvent 

Triarylamine 

978 

Acetone 

Tetraaryldiamine 

949 

Acetone 

Metal  Dithiolene 

836 

Methylene  Chloride 

Indolium  Iodide 

775 

Methanol 

Light 

Brown 

Desert 

Tan 

Urban 

Tan 

Dark 

Brown 


Figure  1.  50/50  NYCO  Ripstop  fabric  considered  for  dye  application.  Fabric  made  of  50% 
Nylon  /  50%  Cotton  blend. 
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Figure  2.  Experimentally  measured  reflectance  spectra  for  fabric  samples  containing  various 
concentrations  of  triarylamine  dye  as  a  function  of  wavelength.  A  control  spectrum  for  fabric 
without  dye  is  shown  for  comparison. 


X  (nm) 


Figure  3.  Experimentally  measured  reflectance  spectra  for  fabric  samples  containing  various 
concentrations  of  tetraaryldiamine  dye  as  a  function  of  wavelength.  A  control  spectrum  for  fabric 
without  dye  is  shown  for  comparison. 
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Figure  4.  Experimentally  measured  reflectance  spectra  for  fabric  samples  containing  various 
concentrations  of  transition  metal  dithiolene  dye  as  a  function  of  wavelength.  A  control  spectrum 
for  fabric  without  dye  is  shown  for  comparison. 


Figure  5.  Experimentally  measured  reflectance  spectra  for  fabric  samples  containing  various 
concentrations  of  indolium  iodide  dye  as  a  function  of  wavelength.  A  control  spectrum  for  fabric 
without  dye  is  shown  for  comparison. 
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(6a) 


(6b) 


(6c) 
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(6d) 


(6e) 

Figure  6.  Modeled  versus  measured  reflectance  spectra  for  fabrics  containing 
concentrations  of  triarylamine  dye  as  a  function  of  wavelength. 
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Figure  7.  Modeled  versus  measured  reflectance  spectra  for  fabrics  containing 
concentrations  of  tetraaryldiamine  dye  as  a  function  of  wavelength. 
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Figure  8.  Modeled  versus  measured  reflectance  spectra  for  fabrics  containing 
concentrations  of  transition  metal  dithiolene  dye  as  a  function  of  wavelength. 
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Figure  9.  Modeled  versus  measured  reflectance  spectra  for  fabrics  containing 
concentrations  of  indolium  iodide  dye  as  a  function  of  wavelength. 
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APPENDIX 


Dielectric  Response  Functions  for  Dyes  Obtained  by  Inverse  Analysis 


Supplementary  information  is  given  here  for  the  inverse  analyses  described  in  references  12 
and  13,  which  are  for  determination  of  estimated  dielectric  response  properties  of  triarylamine 
(Dye  978),  tetraaryldiamine  (Dye  949),  transition  metal  dithiolene  (Dye  836)  and  indolium 
iodide  (Dye  775)  [2],  For  transmission  through  a  sample  of  finite  thickness,  a  consistent 
parametric  representation  of  the  transmissivity  function  is  given  by 
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where  fri,  fa  and  Cri,  Cu  represent  an  optimal  set  of  basis  functions  and  their  weighting 
coefficients,  respectively,  for  a  given  set  of  transmission-spectra  measurements.  The  basis 
functions  describe  the  permittivity  function  within  the  region  bounded  by  wavelengths  Ami-„  and 
kmax,  respectively.  The  wavelengths  Xn  and  quantities  Tc  specify  constraint  conditions  on  the 
transmissivity.  For  parametric  representation  of  the  sample’s  dielectric  response,  a  sufficiently 
optimal  set  of  basis  functions  include  Lorentz  oscillators,  as  well  as  other  functions  satisfying  the 
Kramer- Kronig  relations  [7], 
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The  formal  procedure  underlying  the  inverse  method  considered  here  [12,13]  entails  the 
adjustment  of  the  permittivity  function  defined  over  the  entire  wavelength  region.  This  approach 
defines  an  optimization  procedure  where  the  transmission  spanning  the  range  of  wavelengths  is 
adopted  as  the  quantity  to  be  optimized.  The  constraint  conditions  are  imposed  on  the 
transmission  by  minimizing  the  objective  function  defined  by 

N  2 

2t  =  1Mt(K)-tM)  (15) 

1  n= 1 

where  Tm(An)  is  the  measured  or  target  transmission  for  wavelength  A„.  The  quantities  vv„ 
(i n=l,...,N)  are  weight  coefficients  that  specify  relative  levels  of  influence  associated  with 
constraint  conditions  Tm(An).  A  consistent  parametric  representation  (Eq.(14))  of  the  permittivity 
function  (A)  =  esr(A)+iesi( A)  is  that  of  a  linear  combination  of  Lorentz  oscillators,  i.e., 

=  1  +  x„  =  1  ■ +  2  XX K AP  g,„ K, )  (16) 
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and  the  quantities  xsr  an<l  Xs,  are  the  real  and  imaginary  parts  of  the  susceptibility  function.  The 

procedure  for  inverse  analysis  defined  by  Eqs.  (7)  through  (18)  entails  adjustment  of  the 
parameters  Aj,  gj  and  Aj  defined  over  the  entire  wavelength  range  of  interest.  This  procedure,  by 
inclusion  of  Eq.(15),  is  based  on  least-squares  parameter  optimization. 
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Figure  10.  (a)  Transmission  spectra  as  a  function  of  triarylamine  dye  concentration  in  acetone, 
where  each  spectrum  was  both  normalized  and  background-corrected,  (b)  Objective  function  for 
model  parameter  optimization  with  respect  to  measured  transmission  spectra  for  0.03  gpL 
triarylamine  dye  in  acetone,  (c)  Estimated  imaginary  part  of  complex  susceptibility  to  within 
adjustable  scale  factor  for  triarylamine  dye.  (d)  Estimated  real  part  of  complex  susceptibility  to 
within  adjustable  scale  factor  for  triarylamine  dye. 
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Figure  11.  (a)  Transmission  spectra  as  a  function  of  tetraaryldiamine  dye  concentration  in 
acetone,  where  each  spectrum  was  both  normalized  and  background-corrected,  (b)  Objective 
function  for  model  parameter  optimization  with  respect  to  measured  transmission  spectra  for 
0.03  gpL  tetraaryldiamine  dye  in  acetone,  (c)  Estimated  imaginary  part  of  complex  susceptibility 
to  within  adjustable  scale  factor  for  tetraaryldiamine  dye.  (d)  Estimated  real  part  of  complex 
susceptibility  to  within  adjustable  scale  factor  for  tetraaryldiamine  dye. 
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Figure  12.  (a)  Transmission  spectra  as  a  function  of  transition  metal  dithiolene  dye  concentration 
in  methylene  chloride,  where  each  spectrum  was  both  normalized  and  background-corrected,  (b) 
Objective  function  for  model  parameter  optimization  with  respect  to  measured  transmission 
spectra  for  0.03  gpL  transition  metal  dithiolene  dye  in  methylene  chloride,  (c)  Estimated 
imaginary  part  of  complex  susceptibility  to  within  adjustable  scale  factor  for  transition  metal 
dithiolene  dye.  (d)  Estimated  real  part  of  complex  susceptibility  to  within  adjustable  scale  factor 
for  transition  metal  dithiolene  dye. 
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Figure  13.  (a)  Transmission  spectra  as  a  function  of  indiolium  iodide  dye  concentration  in 
methanol,  where  each  spectrum  was  both  normalized  and  background-corrected,  (b)  Objective 
function  for  model  parameter  optimization  with  respect  to  measured  transmission  spectra  for 
0.03  gpL  indiolium  iodide  dye  in  methanol,  (c)  Estimated  imaginary  part  of  complex 
susceptibility  to  within  adjustable  scale  factor  indiolium  iodide  dye.  (d)  Estimated  real  part  of 
complex  susceptibility  to  within  adjustable  scale  factor  for  indiolium  iodide  dye. 
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